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We report a numerical investigation of the Anderson transition in two-dimensional systems with 
spin-orbit coupling. An accurate estimate of the critical exponent v for the divergence of the 
localization length in this universality class has to our knowledge not been reported in the literature. 
Here we analyse the SU(2) model. We find that for this model corrections to scaling due to irrelevant 
scaling variables may be neglected permitting an accurate estimate of the exponent v — 2.73 ±0.02. 



According to the scaling theory of locahsation if 
interactions are neglected, all states are localised in two 
dimensions (2D). Two exceptions predicted by the scal- 
ing theory to this oft recited mantra are the extended 
states which occur at the center of a Landau level in the 
quantum Hall effect (QHE) and the Anderson tran- 
sition which occurs in zero magnetic field if there is a 
significant spin-orbit interaction 

A surprising aspect of the QHE is that the estimate 
v — 2.35 ± .03 1^, ^ for non-interacting electrons of the 
critical exponent v, which describes the divergence of the 
localisation length ^ at the transition, is in close agree- 
ment with the measured value. This is despite the fact 
that interactions are clearly relevant since the dynami- 
cal exponent z is predicted to be two for non-interacting 
electrons while the experimental value is unity. 

Critical phenomena are determined by the symmetry 
of the Hamiltonian and the dimensionality of the system. 
The important symmetries for the Anderson transition 
exhibited by non-interacting electrons are time reversal 
symmetry and spin rotation symmetry. There are three 
universality classes: orthogonal, unitary and symplectic. 
Systems with time reversal symmetry but where spin ro- 
tation symmetry is broken by the spin-orbit interaction 
belong to the symplectic class. 

In this paper we estimate using numerical simulation 
and finite size scaling the exponent v for the Anderson 
transition in the symplectic universality class in 2D. In 
contrast to the QHE, where extended states occur only 
in vanishingly small energy region, the metallic phase ex- 
tends over a finite energy interval. This system is, there- 
fore, a good candidate for the study of 2D quantum phase 
transitions. Early work suggested that metallic phase in 
this model is destroyed when interactions between elec- 
trons are taken into account ||^, , while more recent 
work suggests that this is not so g. 

Recently, there have been numerous reports of the ob- 
servation of a zero magnetic field metal-insulator transi- 
tion in 2D together with measurements of the critical ex- 
ponents which characterize this transition [^. Whether 
or not those experiments indicate the existence of a true 



TABLE I: Published estimates of the critical exponent for 
the 2D symplectic universality class. QIDLL means finite 
size scaling for the quasi-one dimensional localization length, 
ELS energy level statistics and MFSS multifractal finite size 
scaling. The errors quoted are one standard deviation. 
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2.05 ±0.08 
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2.75 ±0.1 


U 
~ 




Evangelou model 


QIDLL 


2.5 ±0.3 
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2.51 ±0.18 


u 




network model 


QIDLL 


2.88 ±0.15 
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Ando model 


MFSS 


2.41 ± 0.24 



metallic phase at zero temperature remains in dispute. 
If there is indeed a transition, the physics driving it and 
whether there is a relation with the transition we study 
here is not yet clear. However, just as was the case for 
the QHE, we believe that an accurate estimate for the 
metal-insulator transition in zero field in non-interacting 
2D systems may prove useful. 

There has been only very limited success in estimating 
the critical exponent with field theoretic methods |l^, 
p7t . In Tabic | we tabulate the estimates of the exponent 
reported in previous numerical studies ^, |l^, |l^, |l^, |l^, 
p^ , p^ . There is considerable variation between these 
estimates. The estimates reported in ^ |l^, |l^, |l^ 
seem to be consistent with a true value of the exponent 
in the range [2.6,2.9]. However, the estimate of |l2| is 
somewhat below this and that of Q is in contradiction 
with the estimates of Q |ljl . 

We use the transfer matrix method ||l^, ^ to estimate 
the localization length A of electrons on an L x quasi- 
ID strip and then extract the critical exponent from a 
finite size scaling analysis of the dependence of A on L and 
disorder. Two important factors limiting the accuracy of 
the estimate of the exponent obtained in this way are 
the accuracy of the data for A and the maximum width 



2 



L for which data are available. The standard error in the 
estimate of A decreases as \J\/Lz, while for a fixed Lz the 
CPU time needed increases as Li^. At the critical point 
A increases linearly with L so that the CPU time needed 
to estimate A to a given accuracy increases as L^. This 
means that it is somewhat easier to improve the accuracy 
of the numerical data than to increase the size of the 
systems simulated. However, as the accuracy of the raw 
data improves, corrections to scaling due to irrelevant 
scaling variables become more important. While such 
corrections can be taken into account pO| ] , the number of 
fitting parameters is increased and correspondingly the 
uncertainty in the estimate of the exponent is increased. 
It is therefore advantageous to choose a model for which 
such corrections are negligible even when the raw data 
are of high accuracy. In this paper we report results for 
an SU(2) model for which this condition is satisfied. 

The Hamiltonian for the SU(2) model describes non- 
interacting electrons on a 2D square lattice with nearest 
neighbour hopping 



H 



E 



V 



E 

<i,j> ,(t,(t' 



(1) 



Here c] ^ (ci.tr) denotes the creation (annihilation) op- 
erator of an electron at the site i with spin a and 
denotes the random potential at site i. We assume a box 
distribution with each uniformly and independently 
distributed on the interval [-W^/2, 1^/2]. The width W 
of the distribution measures the strength of the random- 
ness. The constant V is taken to be the unit of energy, 
V = l. 

The spin-orbit coupling appears in the hopping matrix 
R{i,j) between each pair of nearest neighbours on the 
lattice. These matrices belong to the group SU(2) of 2 x 
2 unitary matrices with determinant one. The hopping 
matrices are parameterised as follows 

j^^^^^^J e»-cosA, e^^-^smP^A ^2) 
\-e "^^''^ sin (3i J e cos/3,;j / 

This matrix describes a rotation of the electron spin in 
three dimensional space. (The Euler angles of this rota- 
tion are related to, but not equal to, the angles a, /3 and 
7.) In the SU(2) model the distribution of these angles is 
chosen so that the R{i,j) are uniformly distributed with 
respect to the group invariant measure (Haar measure) 
on SU(2). This corresponds to a and 7 uniformly dis- 
tributed in the range [0, 27r), and /3 distributed according 
to the probability density. 



Pimp 



sin(2/3)d/3 < /3 < f 
otherwise. 



(3) 



Hopping matrices on different bonds of the lattice are 
statistically independent. Periodic boundary conditions 



are imposed in the transverse direction. The necessary 
calculations are carried out using quaternion arithmetic 
[Q which halves the required number of multiplications 
compared with an implementation using complex arith- 
metic. 

Some of the physics of the SU(2) model can be under- 
stood by comparing it with the Ando model which has 
been adopted in [|[ ^ ^ In the Ando model as 
the electron propagates through the material its spin pre- 
cesses at a rate and about an axis which depend on the 
electrons wave number. When scattered by the random 
potential, the rate and the axis about which the electron's 
spin rotates changes. This leads to a diffusive motion of 
the spin with an associated spin relaxation length. Quan- 
tum interference between time reversed electron trajec- 
tories longer than this length produces the weak anti- 
localization effect |2^. Motivated by the conjecture that 
the spin relaxation length might be an important irrele- 
vant length scale, we adopted the SU(2) model where the 
uniform distribution of the hopping matrices on SU(2) 
ensures that the spin relaxation length is the shortest 
possible. Doing so we do indeed find that corrections 
due to irrelevant scaling variables can be neglected. 

To determine the critical exponent v, critical disorder 
Wc and other critical properties of the transition we fit 
dependence of A = X/L on the system size L and the 
disorder W, or when = on the Fermi energy Ep, to 
a one parameter scaling law of the form 



In A = F{iP)^ 



(4) 



Here is the relevant scaling variable. We expand the 
scaling function as a power series 



F (x) = lnAc + X + a2x'^ + . . . + a„„a;'^ 



(5) 



terminating the expansions at order no- To allow for 
non-linearity of the scaling variable, the scaling variable 
is approximated by an expansion in terms of the dimen- 
sionless disorder w — (Wc — W)/Wc, where Wc is the 
critical disorder separating the insulating and metallic 
phases. (If = we set w = (Ep - Ec)/Ec.) The 
growth of the relevant scaling variable with system size 
is described by the critical exponent v 



(6) 



where we terminate the expansion at order n^. This 
same exponent describes the divergence of the localiza- 
tion (correlation) length 



(7) 



where we terminate the expansion at order n^. The ab- 
solute scales ^± of the localization length on either side 
of the transition are not determined in this analysis, so 
we set them both to unity for simplicity. The linear co- 
efficient in the expansions of F is set to unity, as shown, 
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to eliminate some redundancy in the definition of the 
fitting parameters. The total number of parameters is 
Np = no + + 2. 

The best fit is determined by minimizing the statis- 
tic. The quality of the fit is assessed with the goodness 
of fit probability Q. Confidence intervals for the fitted 
parameters are estimated using a Monte Carlo method 
p2| . This involves using the model and the best esti- 
mates of the fitting parameters to generate an ideal data 
set. From this data set a large ensemble of synthetic data 
sets is generated by adding random errors, with a vari- 
ance equal to that of the error of the corresponding data 
point, to the ideal data set. Fitting of the ensemble of 
synthetic data sets produces a distribution for the critical 
parameters from which confidence intervals and the good- 
ness of fit are estimated. This procedure is standard and 
systematic but does not take into account any unknown 
systematic effects that might only be discernible for very 
much larger systems. Of course, this caveat applies to 
almost any numerical estimate of a critical exponent. 

Before turning to the estimate of the critical exponent 
we sketch the phase diagram in Figure 0. The fi gure IS 
based on data for systems with sizes ranging from L — 
8 to L = 32 for a number of Fermi energies between 
Ep = and Ep = 3 for which the critical disorder Wc 
was estimated, and on data for = for which the 
critical energy Ec was estimated. In the absence of a 
random potential i.e. when W = 0, the Hamiltonian may 
have chiral symmetry and, in addition to the transition at 
finite Ec, a critical state may also be present at the band 
center Whether this possibility is realised depends 
on the boundary conditions and whether the number of 



TABLE II: The details of various fits to the numerical data 
for the SU(2) model. The fit is to Nd = 230 data points in 
the range 5.2 < W < 6.7 and 0.2 < InA < 0.9. 



no Np Q Wc In Ac u 
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0.2 


5.952± 


.002 


0.612± 


.001 


2.74± .01 
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0.4 


5.952± 


.002 


0.612± 


.001 


2.73± .02 
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0.3 


5.952± 


.002 


0.612± 


.001 


2.73± .02 
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0.3 


5.952± 


.002 


0.612± 


.001 


2.74± .03 
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3 


9 


0.3 


5.952± 


.002 


0.612± 


.001 


2.73± .03 



sites is even or odd. Since chiral symmetry is broken by 
a random potential it does not affect our estimation of v 
below. 

To estimate the critical exponent accurately more ex- 
tensive simulations were carried out for a single energy 
Ep — 1. The numerical data are presented in Figures ^ 
and ||. Data with an accuracy of 0.1% are available for 



TABLE III: The variation of the estimates of the critical pa- 
rameters for the SU(2) model as data for smaller systems sizes 
are progressively excluded from consideration. Here no = 3 
and n^p = 2. 

Nd Q Wc InAe ly 

8 < L < 96 230 0.4 5.952± .002 0.612± .001 2.73± .02 
16 < i < 96 169 0.6 5.953± .003 0.611± .002 2.75± .02 
32 < L < 96 113 0.5 5.954± .005 0.611± .003 2.71± .04 
64 < L < 96 64 0.8 5.96± .02 0.60± .02 2.8± .2 
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TABLE IV: The variation of the estimates of the critical 
parameters for the SU(2) model as the range of disorder under 
consideration is progressively narrowed. 



W Na no 


Q 




In Ac 








[5.2,6.7] 230 3 


2 0.4 


5.952± .002 


0.612± 


.001 


2.73± 


.02 


[5.5,6.4] 175 2 


2 0.7 


5.953± .003 


0.612± 


.001 


2.75± 


.02 


[5.8,6.1] 65 1 


1 0.5 


5.950± .002 


0.613± 


.002 


2.72± 


.07 



system sizes L = 8, 16 and 32, with accuracy 0.2% for 
L = 64, and 0.4% for L = 96. This required Lz of the or- 
der of 10^ to 10^ depending on the size, the disorder and 
the accuracy. When fitting the data the intervals of W 
and In A to consider must be decided. The exact choice 
is not particularly important provided all data are in the 
critical regime. A larger interval of In A requires a higher 
order of expansion in Eq. ^, while a larger interval of W 
requires a higher order expansion in Eq. |6[ The results of 
the finite size analysis are presented in Tables |l[ [II and 
IV. A number of fits of the numerical data are possible 
but, as can be seen by referring to Table ||, all yield con- 
sistent results. The estimates of the critical parameters 
are also stable against restriction of the system sizes un- 



der consideration, see Table III , an d ag ainst a narrowing 
of the range of disorder, see Table [V, The lines shown 
in Figure y| and ^ correspond to a fit with ng = 3 and 

— 2. To demonstrate single parameter scaling graph- 
ically we re-plot the data as a function of L/^ in Figure 
^; metallic and insulating branches are clearly visible. 
These are described by related scaling functions and 

where 



F± (x) ^ In Ac±a; 



-a2X 



2/v 



(8) 



In summary, we have studied the Anderson transition 
in the 2D symplectic universality class and estimated the 
critical exponent v. We find ly = 2.73 ± .02, where the 
error is a 95% confidence interval. Our result is consis- 
tent with the estimates reported in |ll|, |l^ 
but not with those of Analyses based on energy 

level statistics, such as [p_2[, also seem to have a tendency 
to produce estimates which are lower than those of the 
transfer matrix method for the 3D orthogonal universal- 
ity class . On the other hand, in our opinion, the error 
bar for the exponent claimed in [B is too optimistic. 



* slevin@phys.sci.osaka-u.ac.jp 
[1] E. Abrahams, P. W. Anderson, D. C. Licciardello, and 

T. V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 
[2] B. Huckestein, Rev. Mod. Phys. 67, 357 (1995). 
[3] S. Hikami, A. I. Larkin, and Y. Nagaoka, Prog. Theor. 

Phys. 63, 707 (1980). 
[4] T. Ando, Phys. Rev. B 40, 5325 (1989). 




TTTn| 1 I I lllll| 1 I llllll| 1 I llllll| 1 I llllll| 1 I llllll| 1 I 

10"' 10"' 10"' 10"' 10"' 10"' 10"' 10"' 10"' 



FIG. 3: When plotted versus L/^ the data divide naturally 
into an upper metallic branch F+ (i/^) and a lower insulating 
branch F- (L/C). 



[5 

[6: 
[7: 

[8 
[9 

[lo; 



[11 

[12; 

[13 
[14 

[is: 
[16: 

[17 

[is: 

[19 
[20 

[21 
[22: 



[23: 

[24: 



B. Huckestein and B. Kramer, Phys. Rev. Lett. 64, 1437 
(1990). 

B. L. Altshuler and A. G. Aronov, Solid State Commun. 
46, 429 (1983). 

C. Castellani, C. D. Castro, G. Forgacs, and S. Sorella, 
Solid State Commun. 52, 261 (1984). 

C. Castellani, C. D. Castro, and P. A. Lee, Phys. Rev. B 

57, R9381 (1998). 

E. Abrahams, S. V. Kravchenko, and M. P. Sarachik, 
Rev. Mod. Phys. 73, 251 (2001). 

U. Fastenrath, G. Adams, R. Bundschuh, T. Hermes, 
B. Raab, I. Schlosser, T. Wehner, and T. Wichmann, 
Physica A 172, 302 (1991). 

S. N. Evangelou, Phys. Rev. Lett. 75, 2550 (1995). 

L. Schweitzer and I. K. Zharekeshev, J. Phys. Cond. 

Matt. 9, L441 (1997). 

R. Merkt, M. Janssen, and B. Huckestein, Phys. Rev. B 

58, 4394 (1998). 

K. Minakuchi, Phys. Rev. B 58, 9627 (1998). 

K. Yakubo and M. Ono, Phys. Rev. B 58, 9767 (1998). 

S. Hikami, Prog. Theor. Phys. Suppl. 107, 213 (1992). 

E. Brezin and S. Hikami, Phys. Rev. B 55, R10169 

(1997). 

J.-L. Pichard and G. Sarma, J. Phys. C 14, L127 (1981). 
A. MacKinnon and B. Kramer, Z. Phys. B 53, 1 (1983). 
K. Slevin and T. Ohtsuki, Phys. Rev. Lett. 82, 382 
(1999). 

G. Bergmann, Solid State Commun. 42, 815 (1982). 
W. H. Press, A. A. Teukolsky, W. T. Vetterling, and B. P. 
Flannery, Numerical Recipes in Fortran (Cambridge Uni- 
versity Press, Cambridge, 1992). 

P. Brouwer, C. Mudry, and A. Furusaki, Phys. Rev. Lett. 
84, 2913 (2000). 

M. L. Ndawana, R. A. Romer, and M. Schreiber, Eur. 
Phys. J. B 27, 399 (2002). 



